Layer— and bulk roton excitations of 4 He in porous media 
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We examine the energetics of bulk and layer-roton excitations of 4 He in various porous medial such 
as aerogel, Geltech, or Vycor, in order to find out what conclusions can be drawn from experiments 
on the energetics about the physisorption mechanism. The energy of the layer-roton minimum 
depends sensitively on the substrate strength, thus providing a mechanism for a direct measurement 
of this quantity. On the other hand, bulk-like roton excitations are largely independent of the 
interaction between the medium and the helium atoms, but the dependence of their energy on the 
degree of filling reflects the internal structure of the matrix and can reveal features of 4 He at negative 
pressures. While bulk-like rotons are very similar to their true bulk counterparts, the layer modes 
are not in close relation to two-dimensional rotons and should be regarded as a third, completely 
independent kind of excitation. 



INTRODUCTION 



Collective excitations of superfluid helium confined in 
various porous media have been studied by neutron scat- 
tering since early 90's, and by now a wealth of informa- 
tion about helium in aerogel, Vycor and Geltech has been 
collected 0-0. Aerogel is an open gel structure formed 
by silica strands (S1O2). Typical pore sizes range from 
few A to few hundred A , without any characteristic pore 
size. Vycor is a porous glass, where pores form channels 
of about 70 A diameter. Geltech resembles aerogel, ex- 
cept that the nominal pore size is 25 A ||J. 

Liquid 4 He is adsorbed in these matrices in the form 
of atomic layers, the first layer is expected to be solid; 
on a more strongly binding substrate, such as graphite, 
one expects two solid layers. Energies and lifetimes of 
phonon-roton excitations for confined He are nearly 
equal to their bulk superfluid 4 He values but dif- 
ferences appear at partial fillings. The appearance of 
ripplons is tied to the existence of a free liquid surface; 
neutron scattering experiments show clearly their pres- 
ence in adsorbed films Elll III with few layers of 
helium. 

An exclusive feature of adsorbed films is the appear- 
ance of "layer modes" . The existence of such excitations 
has been proposed in the seventies [TiLll"^ from theoreti- 
cal calculations of the excitations of two-dimensional 4 He 
and comparison with specific heat data. Direct experi- 
mental evidence for the existence of collective excitations 
below the roton minimum has first been presented by 
Lauter and collaborators [l3L Hf^ | , identification of these 
excitations with longitudinally polarized phonons that 
propagate in the liquid layer adjacent to the substrate 
has been provide d by microscopic calculations of the ex- 
citations of films 13 ■ 

In an experimental situation, the topology gives rise 
to non-uniform filling of the pores. But from the the- 
oretical point of view different materials are character- 
ized solely by their substrate potentials, because as long 
as the wavelength of the excitation in concern is much 
shorter than any porosity length-scale, the topology of 
the confining matrix is immaterial. We therefore exam- 
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FIG. 1: The density profiles of the first three layers are shown 
as a function of the depth D of the substrate potential. The 
substrate is at z < 0. 



ine the energetics of the layer-roton as a function of the 
substrate-potential strength which determines, in turn, 
the areal density in the first liquid layer. For that pur- 
pose, we have carried out a number of calculations of 
the structure of helium films as a function of potential 
strength. The microscopic theory behind these calcula- 
tions is described in Ref. O ur model assumes the 
usual 3-9 potential 



4gf 
27D 2 



z 3 



(1) 



we have varied the potential strength D from 8 K to 50 K 
and the range C 3 from 1000 K A 3 to 2500 KA 3 . In all 
cases, we have considered rather thick films of an areal 
density of 0.45 A~ 2 . Fig. ^ shows density profiles for 
these potential strengths close to the substrate; the den- 
sity profiles are practically independent of the potential 
range C 3 . 
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II. THEORY OF EXCITATIONS 

To introduce excitations to the system one applies a 
small, time-dependent perturbation that momentarily 
drives the quantum liquid out of its ground state. Gener- 
alizing the Feynman-Cohen wave function |20j | , we write 
the excited state in the form 



!*(*)> = 



t/h e ^SU(t) 

|(* |e ^( 4 )|* )] 1/2 



(2) 



where \^o) is the exact or an optimized variational 
ground state, and the excitation operator is 



(3) 



The time-dependent excitation functions 
5u n (ri, . . . , r n ; t) are determined by an action prin- 
ciple 



H 



d 

ih— + U cxt {t) 



*(t))=0, (4) 



where U ex t(t) is the weak external potential driving the 
excitations. The truncation of the sequence of fluctuat- 
ing correlations Su n in Eq. defines the level of ap- 
proximation in which we treat the excitations. One re- 
covers the Feynman theory of excitations [2lJ for non- 
uniform systems by setting <5m 2 (i"i, • ..r n ;t) — for 
n > 2. The two-body term (5u 2 (ri, r 2 ; t) describes the 
time-dependence of the short-ranged correlations. It is 
plausible that this term is relevant when the wavelength 
of an excitation becomes comparable to the interparticle 
distance. Consequently, the excitation spectrum can be 
quite well understood [23l Y2A l25| by retaining only the 
time-dependent one- and two-body terms in the exci- 
tation operator @. The simplest non-trivial implemen- 
tation of the theory leads to a density-density response 
function of the form 0] 

x(rY,u/)= (5) 

St 

where the (f>^ (r) are Feynman excitation functions, and 



Gst(w) = [h{u - uj s + ie]S s t + Sst(w)]" 



(6) 



the phonon propagator. The fluctuating pair correlations 
give rise to the dynamic self energy correction |l7j |. 



V (s)y(t) 
y mn * mn 



2 ^ h(uj ni + UJ n - uj) 



(7) 



Here, the summation is over the Feynman states to, n; 
they form a partly discrete, partly continuous set due 
to the inhomogeneity of the liquid. The expression for 



(s) 

the three-phonon coupling amplitudes Vmn can be found 
in Ref. 17. This self energy renormalizes the Feynman 
"phonon" energies fiuj n , and adds a finite lifetime to 
states that can decay. The form of the self energy given in 
Eq. is the generalization of the correlated basis func- 
tions (CBF) pH. |24| theory to inhomogeneous systems. 
As a final refinement of the theory, we scale the Feynman 
energies w„ appearing in the energy denominator of the 
self energy given in Eq. Q such that the roton minimum 
of the spectrum used in the energy denominator of Eq. 
(J7J agrees roughly with the roton minimum predicted by 
the calculated SQc,u>). This is a computationally simple 
way of adding the self energy correction to the excitation 
energies in the denominator of Eq. Q . We shall use this 
approximation for the numerical parts of this paper. 



III. LAYER EXCITATIONS 

Layer phonons are identified by a transition density 
that is localized in the first liquid layer of the system. 
They appear in the dynamic structure function 5(k, uS) 
as a peak below the roton minimum. A grayscale map of 
a typical dynamic structure function is shown in Fig. |2 
we have for clarity chosen a momentum transfer parallel 
to the substrate; neutron scattering at other angles would 
broaden the roton minima [Isj |. The figure shows in fact 
one bulk and two layer-roton minima, but the higher 
one, which corresponds to an excitation propagating in 
the second liquid layer, has an energy too close to the 
bulk roton to be experimentally distinguishable. 

The transition densities corresponding to the three 
pronounced excitations at k = 1.8 A -1 are depicted in 
Fig. |3 Clearly, the two "layer-modes" are actually lo- 
cated in the two first layers adjacent to the substrate 
whereas the "bulk" mode is spread throughout the sys- 
tem. However, the figure also shows that the notion that 
the wave propagates in the first or the second layer is 
also not quite accurate: The lowest mode also has some 
overlap with the second layer, but especially the second 
mode spreads over both layers. 

We have carried out two independent calculations of 
the two-dimensional roton excitation: First, we calcu- 
lated the roton energy as a function of the density for a 
rigorously two dimensional liquid. We can assess the ac- 
curacy of our predictions with the shadow-wave-function 
calculation of 5(k, uj) of Ref. |2y, who obtained a ro- 
ton energy of 5.67 ± 0.2 K at the equilibrium density of 
n = 0.0421 A -2 . Second, we have calculated the dynamic 
structure function S(k., to) in the relevant momentum re- 
gion for the above family of substrate potentials. The 
results are compiled in Fig. 0] where we also collect sev- 
eral experimental values. 

Although exactly the same method has been used for 
the computation of the purely two-dimensional system 
and for the films, the results are quite different. We have 
obtained for the film calculation an effective layer den- 
sity by integrating the three-dimensional densities shown 
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FIG. 2: The figure shows the map of the dynamic structure 
function S(k, lu) for a 4 He film for the potential strength D = 
32 K. The two layer-rotons, the bulk roton, and the ripplon 
are indicated by arrows. 
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FIG. 3: The transition densities of the three lowest excitation 
are shown for D = 32 K and C3 = 1500 KA 3 , normalized to 
the same maximum value. For comparison, the density profile 
of the film is shown as a gray-shaded area. 
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FIG. 4: The figure shows the energy of the roton in two di- 
mensions, (long-dashed line) and the energy of the layer roton 
in a film (solid line) as a function of areal density. The areal 
densities of the layers were obtained by integrating the den- 
sity up to the first minimum, the upper horizontal axis shows 
the corresponding values of the well-depth D. Also shown is 
the energy of a two-dimensional roton obtained with shadow 
wave functions |2(| at the density of 0.0421A" 2 . The short- 
dashed horizontal lines show experimental values of the roton 
energy on aerogel |l[lQ] (7.39 and 7.14 K), Vycor0G| (7.23 
K), and Geltech|jJ (6.71 K), their energies are marked on the 
right margin. 

111 Fig. □ to the first minimum. This is evidently not 
very well defined for the weakly bound systems, but it 
is not legitimate cither for the case of strong binding 
where the first layer is well defined. In fact, the integrated 
density for the strongest substrate is 0.08 A -2 , which is 
well beyond the solidification density of the purely two- 
dimensional system. Evidently, the zero-point motion in 
z direction can effectively suppress the phase transition. 
We make therefore three conclusions: (i) The position of 
the layer roton minimum is indeed a sensitive measure for 
the strength of the substrate potential, (ii) purely two- 
dimensional models are manifestly inadequate for their 
understanding, and, hence, (iii) purely two-dimensional 
models are also questionable for interpreting thermody- 
namic data of adsorbed films. 



IV. BULK EXCITATIONS 

With one exception the bulk roton energy in porous 
media have been reported to be practically identical to 
that in the bulk liquid, Ref. reports a slight increase 
of the roton energy in aerogel at partial filling. A roton 
energy above the one of the bulk liquid can be explained 
by assuming that the density of the helium liquid in the 
medium is below that of the bulk liquid. This can, in 
turn, be qualitatively explained by the cost in energy to 
form a surface. 

To be quantitative, we have performed calculations of 
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FIG. 5: The figure shows the energy of the bulk roton in 
a gap of 25 A width as a function of the two-dimensional 
coverage and the corresponding average density. The points 
of zero pressure, the low-density and the high-density spin- 
odal points, and the bulk chemical potential are indicated as 
points. The dashed line shows the roton energy of a purely 
3D calculation. 

the energetics and structure of 4 He in a gap between at- 



tractive silica walls [18j and obtained the energy of the 
bulk roton (c.f. Fig. [21 as a function of filling. Fig. 
shows, as a typical example, the roton energetics in 
in a gap of 25 A width. The independent parameter is 
the areal density n, the corresponding three-dimensional 
density was obtained by averaging the density profile over 
the full volume. It is seen that the equilibrium density 
is well below the bulk value. In other words, the roton 
energy in a confined liquid should correspond to the one 
of a liquid that would, without confinement, have a nega- 
tive pressure. The energy increase of the roton minimum 
found in this model is about 0.5 K, which is consistent 
with the experiments of Ref. ||J 

To verify this interpretation of the data, it would be 
very useful to have comparable measurements for porous 
media with a more uniform distribution of pore sizes. 
In particular, comparably small pores should allow to 
densities that are even below the bulk spinodal density 
|27| , thus facilitating experiments on 4 He in density areas 
that were up to now inaccessible. 
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